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ABSTRACT 


V 


An  experiment  was  conducted  in  the  Florida  Keys  from  2  -  19  July  1990  to  test 
the  performance  of  the  Naval  Oceanic  Vertical  Aerosol  Model  (NOVAM)  in  a  weak- 
convective  regime.  Meteorological  data  collected  by  aircraft  and  boat  was  used  to  gen¬ 
erate  the  surface  and  vertical  profile  information  files  required  by  NOVAM,  Using  this 
information,  NOVAM  predicts  the  aerosol  extinction  (/c«r')  for  a  vertical  cross-section 
of  the  atmosphere.  Aircraft-observed  aerosol  extinction  profiles  were  also  obtained. 
Comparisons  between  observed  and  NOVAM  aerosol  extinction  profiles  revealed  major 
deviations  above  the  cloud  top.  From  the  surface  to  the  top  of  the  cloud  layer,  NOVA.M 
generally  did  an  excellent  job  in  predicting  profile  shape,  with  the  magnitude  of  aerosol 
extinction  tied  to  the  extinction  matched  at  the  surface.  In  a  few  cases,  observed  ex¬ 
tinction  i;;. Teased  more  rapidly  than  NOVAM  predicted  extinction  from  the  surface  to 
the  base  of  the  cloud  layer.  This  is  attributed  to  rain  scavenging  associated  with 
thunderstoim  activity  in  tne  area.  Comparison  between  different  aerosol  extinction 
profiles  revealed  much  spatial  and  temporal  variation  that  was  verified  by  Lidar  profiles 
of  atmospheric  structure.  Thunderstorm  activity,  multiple  cloud-layers,  and  the  spatial 
variation  in  the  atmospheric  structure  have  led  to  a  hypo.hesis  that  deep-convection  was 
responsible  for  the  m.ajor  differences  between  observed  and  predicted  aerosol  e,\tinction 
profiles.  If  this  is  the  case,  a  simple  modificaiion  to  the  weak-convective  model  may  be 
made  to  apply  to  a  deep-convection  model. 
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I.  INTRODUCTION 


Successful  employment  of  modem  weapon  systems  require  the  continuous  ability  to 
^  use  wavelengths  in  the  electro-magnetic  (EM)  spectmm.  This  is  the  fundamental  con¬ 

cept  of  electronic  warfare  and  such  importance  has  been  placed  on  this  mission  area  that 
,  the  following  maxim  was  developed;  The  side  that  controls  the  electro-magnetic  (EM) 

spectrum  controls  the  outcome  of  any  conflict  in  modem  war  or  global  politics'.  Oper¬ 
ation  Desert  Storm  was  successful  for  many  reasons,  one  of  which  was  the  allied  forces 
superiority  in  the  maintained  use  of  systems  operating  at  various  wavelengths  through¬ 
out  the  EM  spectrum.  Even  with  this  superiority,  there  is  much  room  for  improvement. 
One  example  is  the  limitation  that  cloud  cover  had  on  operations.  .Another  is  operating 
in  regions  with  extensive  aerosol  content  (such  as  were  caused  by  the  intentional  burning 
of  oil  rigs). 

Many  military  devices  operate  in  the  visible  to  the  far  infrared  (IR)  wavelengths. 
Some  examples  are  forward-looking  infrared  (FLIR),  precision  guided  munitions  and 
other  electro-optical  (EO)  systems.  For  this  reason,  there  is  a  great  deal  of  interest  in 
reliable  estimation  of  EO  instrumentation  performance  for  vertical  and  slant-path  ob¬ 
servations.  For  systems  such  as  those  above,  this  requires  knowledge  of  the  vertical 
variation  of  aerosol  scattering  and  absorption  in  the  wavelengths  of  operation. 

LOWTRAN  7  (Low  resolution  transmittance  model  version  7)  calculates  atmo¬ 
spheric  scattering  and  absorption  caused  by  aerosols  and  molecules  along  a  non- 
hemogeneous  path.  Because  there  isn't  one  aerosol  model  that  predicts  scattering  and 
absorption  the  best  under  all  conditions,  there  are  diflerent  aerosol  models  a  user  of 
LOWTR/\N  may  choose.  For  maritime  applications,  the  present  model  is  the  Navy 
.Aerosol  Model  (N.A.M).  NAM  attempts  to  relate  aerosol  size  distributions  to 
meteorological  parameters  such  as  relative  humidity,  wind  speed,  and  visibility,  and  then 
«.  calculate  the  optical  properties  of  the  modeled  aerosol.  Most  empirical  models,  such  as 

N.AM,  do  not  predict  the  vertical  distribution  of  aerosol  required  to  extend  extinction 
•  predictions  to  higher  levels.  For  the  maritime  regime,  this  led  to  the  development  of  the 

Naval  Oceanic  Vertical  .\erosol  .Model  (NOV.A.M). 

The  purpose  of  NOV.AM  (Gathman,  1989)  is  to  determine  the  vertical  distribution 
of  the  aerosol  size  spectra  in  the  marine  environment,  used  to  obtain  the  associated  op¬ 
tical  and  infrared  properties  along  slant  paths  within  the  marine  boundary  layer.  The 
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means  used  to  achieve  this  goal  are  to  integrate  various  marine  aerosol  studies  of  a  di¬ 
verse  nature  into  a  single  model.  This  integration  has  led  to  the  develop  nc-nt  of  sub¬ 
models  within  NOVA.M,  one  of  which  is  selected  based  on  ttic  parameters  vHst  3rc  input* 
This  thesis  will  investigate  the  performance  of  XOVAM  under  corsditions  of  weak 
cumulus  convection.  These  conditions  existed  during  the  KEY90  experiment  which  took 
place  from  2  July  1990  through  19  July  1990.  near  Marathon,  Florida  on  the  Florida 


II.  BACKGROUND 


A.  EXTINCTION  PARAMETERS 

In  this  section,  fundamental  concepts  and  terms  that  form  the  basi:  of  the 
LOWTRAN  model,  and  led  to  the  development  of  NAM  and  NOVAM  are  presented. 
Lambert,  Beer  and  Bouguet  established  a  relationship  between  attenuated  radiation  and 
the  incident  radiation  for  a  homogeneous  extinction  as  follows: 

I=I,itxp{-bx)  (1) 


where  I  is  the  attenuated  radiation, 
is  the  incident  radiation, 
b  is  the  total  extinction  coefllcient, 

X  is  the  path  length. 

To  determine  the  total  e,xtinction  coefficient,  b,  knowledge  that  extinction  is  caused  by 
absorpfion  and  scattering  by  molecules  and  aerosol  (dry  or  liquid)  contributing  linearly 
is  applied  as  follows: 


^  ~  ^ma  "b  ^ms  "b  ^aa  "b  ^as  (^) 

w'here  b„,  =  extinction  due  to  molecular  absorption, 

b„,  =  extinction  due  to  molecular  (Rayleigh)  scattering, 
b„  =  extinction  due  to  aerosol  absorption, 
b„  =  extinction  due  to  aerosol  (Mie)  scattering. 

In  the  visible  region  of  the  spectrum,  b„,  and  b„,  are  the  dominant  extinction  com¬ 
ponents.  In  the  IR  regions,  b„,  and  b„  are  the  primary  components  causing  lowered 
transmittance  (the  ratio  of  I  to  /„).  Figure  2  on  page  5  shows  the  effect  molecular  ab¬ 
sorption  has  on  transmittance  over  a  horizontal  1.82  km  path  for  wavelengths  up  to  15 
^im.  This  figure  clearly  demonstrates  that  atmospheric  total  extinction  is  wavelength 
dependant  due  to  at  least  its  b„,  component. 

Equation  I  assumes  a  homogeneous  total  extinction  throughout  the  atmosphere. 
Atmospheric  studies  have  shown  this  assumption  to  generally  be  valid  for  horizontal 
paths  (such  as  Figure  2),  however,  a  model  needed  to  be  developed  to  predict  the  vari- 
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Figure  2.  Atmospheric  Transmittance  (0  -  15  microns):  Atmospheric 

Transmittance  measured  over  a  1.82  km  horizontal  path  at  sea  level 
(Hudson,  1969). 

ations  of  extinction  in  the  vertical  required  to  perform  slant-path  calculations. 
LOWTRAN  is  the  model  that  resulted  from  this  requirement.  LOWTRAN  was  de¬ 
signed  to  calculate  the  total  extinction  as  a  function  of  altitude,  and  produce  as  output 
the  radiance  and/or  the  transmittance  for  a  specified  path  through  the  atmosphere. 

B.  THE  LOW  RESOLUTION  TRANSMITFANCE  MODEL  (LOWTRAN) 
LOWTRAN's  radiance  calculations  account  for  contributions  from  atmospheric 

self-emission,  solar  and/or  lunar  radiance  single  scattered  into  the  path,  direct  solar 
irradiance  through  a  slant  path  to  space  and  multiple  scattered  solar  and/or  sclf- 
emmission  radiance  into  a  path.  LOWTRAN  7  (the  latest  version)  incorporates  separate 
molecular  profiles  for  all  major,  as  well  as  13  minor  and  trace  gases  found  in  the  at¬ 
mosphere.  In  addition,  various  aerosol,  cloud  and  rain  models  are  included  to  allow  the 
user  to  select  the  models  that  most  closely  represent  the  region  of  interest.  The  user  also 
has  the  option  to  select  from  six  reference  atmospheric  profiles,  or  define  a  new  atmo¬ 
spheric  profile.  A  more  extensive  description  of  LOWTRAN  may  be  found  in  ONTAR 
corp.  (1990).  (Kneizys  et  al.,  1988) 

C.  THE  NAVY  AEROSOL  MODEL  (NAM) 

NAM  represents  one  of  the  several  empirically  derived  models  contained  in 
LOWTRAN  and  was  designed  to  estimate  the  contribution  of  aerosol  to  the  LO  prop- 
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agating  characteristics  in  the  marine  environment.  This  model  predicts  an  aerosol  size 
distribution  (dX/dr)  within  the  atmospheric  surhice-layer,  nominally  for  10  meters,  for 
specified  wind  speed  (current  and  24-hour  averaged),  visibility,  and  relative  humidity 
(RH)  values.  Three  distinct  types  (modes)  of  aerosols  are  assumed  to  exist  within  the 
size  distribution: 

1.  A  continental  component  contributed  by  a  nearby  land  mass, 

2.  A  stationary  component  affected  by  the  winds  (controlled  by  the  24-hour  averaged 
wind  speed), 

3.  A  fresh  component  caused  by  the  current  wind  over  the  W'ater. 

The  predicted  aerosol  size  distribution  is  obtained  by  summing  the  contribution  of  each 
component  as  follows; 

3 

■f"  =  E  1  ‘"W  -  '"W  -  --Jf )  (3) 

r=l 

where  X  represents  the  number  density  per  unit  radius  of  particles  at  radius  ,r , 

A/  is  the  amplitude  parameter  for  mode  i, 

/  is  the  aerosol  growth  parameter  for  mode  i  (to  be  discussed  in  a  later  section), 
is  the  mode  radius  parameter. 

When  a  standard  relative  humidity  of  80%  is  introduced  into  the  above  equation,  /  is 
replaced  by  the  value  one,  as  will  be  shown  in  the  relative  humidity  discussion.  A  plot 
of  log(dX,  dr)  at  80%  relative  humidity  versus  log(r)  is  shown  in  Figure  3.  In  this  figure, 
Al,  A2,  and  A3  are  the  amplitude-.parameters  as  functions  of  air  mass  type,  24-hour 
averaged  wind  speed,  and  current  wind  speed  respectively;  and  r  is  the  aerosol  particle 
radius  at  the  ambient  relative  humidity.  The  mode  radius  parameter,  r,,,  is  the  radius  that 
the  maximum  amplitude  parameter  value  for  each  mode,  i,  occurs  for  this  80%  relative 
humidity.  (Gathman,  1983) 

XAM's  application  is  to  make  horizontal  path  calculations  of  aerosol  extinction  at 
the  sea  surface  based  on  the  following  relationship  between  dX/dr  and  aerosol  ex¬ 
tinction  (/?„,): 


/^exi  ~  ^aa  "b 


n  dN  2  , 

Qexi  ■  (if  ^ 


(4) 
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Figure  3.  Idealized  three  lognormal  marine  aerosol:  (Gathman,  1983). 

where  Q,„  is  the  Mie  efficiency  coefficient  for  extinction  which  is  a  function  of  wave¬ 
length,  causing  /?„,  to  be  wavelength  dependent. 

To  use  NAM  in  slant  path  calculations,  an  assumed  vertical  size  distribution  must 
be  developed  for  each  mode  based  on  the  size  distribution  NAM  produces  at  10  meters. 
The  present  assumed  size  distribution  profile  has  an  exponential  decrease  with  height. 
NOVAM  is  being  developed  to  provide  a  capability  for  more  accurate  slant-path  aer¬ 
osol  extinction  calculations  in  the  maritime  regime. 

D.  THE  MARITIME  ATMOSPHERIC  BOUNDARY  LAVER  (MABL) 

Vertical  profiles  of  aerosol  depend  on  the  thermodynamic  structure  and  processes 
within  the  MABL.  A  schematic  representation  of  the  MABL  is  shown  in  Figure  4.  In 
this  scheme,  the  cloud  layer  varies  the  most  and,  at  times,  may  not  be  present.  In  cases- 
where  the  large  scale  flow  prevents  cloud  development,  the  MABL  terminates  at  the 
transition  layer.  In  cases  where  very  deep  cloud  development  has  taken  place,  the  top 
of  the  MABL  is  not  well  defined.  Features  such  as  frontogcncsis  and  horizontal 
advection  would  cause  the  one-dimensional  view  of  the  MABL  (as  depicted  in 
Figure  4)  to  be  no  longer  valid. 

For  a  MABL  as  shown  in  Figure  4,  the  following  principles  apply:  (Augstein,  1976) 
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Figure  4.  Schematic  of  dominant  processes  in  a  multi-layered  MABL:  0  =  poten¬ 
tial  temperature.  (Augstein,  1976) 


1.  Radiative  cooling  is  compensated  for  mainly  by  the  adiabatic  warming  of  large- 
scale  subsiding  air  and,  to  a  lesser  extent,  by  the  sensible  heat  input  from  the  ocean. 

2.  Buoyancy  driven  mixing  produces  a  well-mixed  region  above  the  shallow  surface 
layer. 

3.  Kinetic  energy  which  is  generated  by  the  concurrent  heating  and  moistening  of  air 
from  below  is  partly  converted  to  potential  energy  through  entrainment  of  poten¬ 
tially  warmer  air  from  above  through  the  top  of  the  mixed  layer.  This  assumes  that 
buoyant  parcels  overshoot  their  level  of  equilibrium,  thereby  forcing  mass  of  less 
density  downward. 

4.  The  top  of  the  mixed  layer  is  determined  by  a  balance  of  small-scale  mass 
entrainment  from  below  and  the  large-scale  downward  mass  flow. 

5.  Clouds  begin  to  develop  when  convective  plumes  or  bubbles  penetrate  to  their 
condensation  level. 
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III.  THE  NAVAL  OCEANIC  VERTICAL  AEROSOL  MODEL  (NOVAM) 
A.  BACKGROUND 

As  described  previously,  NAM  does  not  predict  the  correct  vertical  distributions 
which  are  required  to  extend  extinction  predictions  to  higher  levels.  NOVAM  is  designed 
to  predict  non-uniform  and  non-logarithmic  extinction  profiles  observed  in  the  MABL 
for  wavelengths  between  .2  and  40  nm.  (Gathman,  1989) 

N’OVAM  is  a  combination  of  empirical  and  physical  models  which  describe  aerosol 
dynamical  behavior.  The  current  version  of  NOVAM  incorporates  the  fr.  jwing  mod¬ 
els: 

1.  N'.A.VI,  which  defines  the  surface  (10  meter)  aerosol  characteristics  (Gathman, 
1983), 

2.  Micro-meteorological  models,  which  describe  the  vertical  distribution  of  aero;  ol  in 
the  well-mi.xed  .MABL  (Fairall  and  Davidson,  1986), 

3.  A  model  that  describes  the  vertical  distribution  of  aerosol  in  the  weak  cum  ilus 
convection  regime  (described  below)  (Davidson  and  Fairall,  1986), 

4.  A  model  describing  the  aerosol  structure  below  stratus  decks  (N'oonkester,  1983), 

3.  .'\  default  boundary  layer  structure  model  of  meteorological  parameters  based  on 
surface  observations  (Gathman,  1989), 

6.  model  that  predicts  aerosol  growih  based  on  relative  humidity  and  the  air  mass 
type  in  the  well-mixed  and  cloud  layers  (Fitzgerald,  1978,  Gerber,  1985). 

Extinction  properties  are  calculated  from  the  generated  aerosol  profiles  using  MIE 
scattering  for  which  two  basic  characteristics  exist: 

1.  The  aerosol  particle  radius  and  the  wavelength  of  incident  energy  are  of  the  same 
order  of  magnitude. 

2.  Scattering  is  generally  concentrated  in  the  forward  direction  (fortt'ard  scatter). 

The  following  particle  source  and  particle  physics  principles  have  been  used  in  the 
development  of  NOVAM.  Large  sea-salt  aerosol  particles  are  produced  from  white 
water  phenomena  at  the  air-sea  interface  (Monahan  et  al.,  1982),  Smaller  aerosol  par¬ 
ticles  are  produced  from  sources  such  as  gas-to-particle  conversion  or  anthropologically 
generated  aerosol.  The  concentration  of  any  particular  aerosol  is  dependant  on  the 
source  strength  of  aerosol  production  and  on  the  mixing  process  as  it  relates  to  scalar 
contaminants.  Regarding  hygroscopic  sea-salt  aerosol,  its  size  is  very  dependant  on  the 
relative  humidity  it  is  immersed  in.  Also,  as  the  hygroscopic  aerosol  picks  up  water  vapor 
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from  the  atmosphere  and  grows  in  size,  it  changes  its  chemical  composition  and  its  index 
of  refra  ction. 

B.  NOVAM  INPUTS 

Tables  1  and  2  show  surface  and  radiosonde  observation  files  (inputs)  required  by 
NOVAM.  The  surface  observed  parameters  are  used  to  perform  the  size  distribution  and 
aerosol  extinction  calculations  at  the  surface  and  for  sub-model  selection,  in  conjunction 
with  some  information  contained  in  Table  2.  The  cloud  type,  present  weather,  and  the 
zonal;  seasonal  category  codes  may  all  be  found  in  separate  tables  as  described  by 
Gathman  (1989).  If  any  of  these  data  items  are  unavailable  for  input,  the  model  calcu¬ 
lates  fiiat  parameter  using  default  information  which  results  in  a  reduced  confidence  in 
aerosol  e.xtinction  predictions. 


Table  I.  SURFACE  OBSERVATION  DATA  FILE 


Position 

.VIeteorological  Data 

1. 

Sea  surface  temperature  (“C) 

2. 

Air  temperature  (°C) 

3. 

Relative  humidity  (RH)  {%) 

4. 

Optical  visibility  (km) 

5. 

Current  real  wind  speed  (m,s) 

6. 

Averaged  wind  speed  (24  h)  (m's) 

7. 

Air  mass  parameter  (see  text) 

8. 

Cloud  cover  (tenths) 

9. 

Cloud  type  [0...9] 

10. 

Surface  IR  ext.  (l/km)  (v'5  I0.6/im) 

11. 

Present  weather  in  standard  code  [0...99] 

12. 

Height  of  lowest  cloud  (m) 

13. 

Zonal'seasonal  category  [1...6] 

.Most  of  the  parameters  listed  in  Table  1  are  self  explanatory',  however,  a  couple 
warrant  further  examination.  The  air  mass  parameter  is  an  indicator  of  the  degree  of 
continental  contaminants  and  may  be  represented  as  follows: 

a.m.p.  =  Rrt/4  +  I 

where  Rn  represents  the  atmospheric  radon  content  (picocuries  per  cubic  cm). 
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Table  2.  CONTENTS  OF  THE  RADIOSONDE  DATA  FILE  PREAMBLE 


Row 

Column  Description 

[!]> 

[1] 

Number  of  actual  radiosonde  observations 

[!]> 

[2] 

Surface  potential  temperature  from  radiosonde 

Cl]‘ 

[3] 

Surface  mixing  ratio  from  radiosonde 

Profile  Characteristics 

[2] 

[1] 

Height  of  base  of  cloud  layer  (CB) 

[2] 

[2] 

Potential  temperature  just  below  CB 

[2] 

[3] 

Mixing  ratio  just  below  CB 

[3] 

[I] 

Units  indicator  for  temperature  (  =  1  if  “Q 

[3] 

[2] 

Potential  temj.ierature  just  above  CB 

[3] 

[3] 

Mixing  ratio  just  above  CB 

[J] 

[1] 

Height  of  the  cloud  layer  top  (CT) 

[•4] 

[2] 

Potential  temperature  just  below  CT 

C'l] 

[3] 

•Mixing  ratio  just  below  CT 

[5] 

[1] 

Units  indicator  for  mixing  ratio  (  =  1  if  g  kg) 

[5] 

[2] 

Potential  Temperature  just  above  CT 

[5] 

ISI 

Mixing  ratio  just  above  CT 

'  Minimum  requirement  for  radiosonde  data  file. 

or  as, 

a.m.p.  =  (9  X  expr  — //4])  +  1 

where  t  equals  the  elapsed  time  (in  days)  it  takes  the  current  air  mass  to  reach  the  point 
of  observation  from  a  distant  land  mass. 

An  air  mass  parameter  of  1  indicates  a  pure  air  mass  (no  contaminants).  When  the  air 
mass  parameter  is  less  than  or  equal  to  5,  the  aerosol  is  assumed  to  be  made  up  of  three 
lognormal  components  (/I,,  /Ij  ,  and  A^)  as  contained  in  the  NAM  description  above. 
For  air  mass  parameters  greater  than  .5,  indicating  the  presence  of  non-soluble  aerosol, 
an  additional  lognormal  component  (Ao)  is  assumed.  This  class  of  non-soluble  aerosol 


is  only  present  close  to  shore  lines  and  has  the  same  mode  radius  as  r,,  shown  in 
Figure  3.  Being  assumed  non-soluble,  the  mode  0  component  of  aerosol  does  not  change 
size  with  changes  in  relative  humidity.  (Gathman.  1939) 

The  remaining  inputs  to  NOVA.M  are  contained  in  a  data  file  describing  the  vertical 
atmosphere,  referred  to  as  the  radiosonde  data  file  in  Gathman,  1989.  Table  2  shows  the 
format  for  the  first  five  rows  that  are  contained  in  this  file,  m.-essary  when  choosing  the 
weak  convective  model.  Figure  5  graphically  illustrates  the  position  in  the  file  by  row 
and  column  designation.  The  remainder  of  the  radiosonde  file  consists  of  radiosonde 
observation  height,  potential  temperature,  and  mixing  ratio  entries.  If  radiosonde  data 
.s  not  available,  a  minimum  of  I  surface  entry  must  be  input  and  a  defaiilt  profile  is  se¬ 
lected. 

After  the  input  parameters  are  translated  by  \OVA'-I  to  predictions  of  mixing, 
source  strengths,  and  size  distributions,  a  prediction  is  provided  of  the  extinction  of  EO 
energy  as  a  function  of  altitude  for  the  wavelength  of  interest.  This  output  is  based  on 
the  dry  aerosol  size  distribution  at  a  particular  height  and  adjusted  for  relative  humidity 
at  that  height. 

C.  SUB-MODELS 

This  section  is  intended  to  give  a  brief  explanation  of  the  different  models  being  de¬ 
veloped  within  NOVA.M.  Figure  6  shows  the  major  decision  points  and  the  track  of 
flow  of  information  in  the  current  version.  There  are  currently  four  sub-models  in 
NOVAM: 

1.  A  mi.xed  boundary  layer  model  (Fairall  &  Davidson,  1986;  Davidson  &  Fairall, 
1986)  where  the  boundary  layer  depth  must  be  less  than  3  km  and  not  stable, 

2.  A  stratus  model  for  winds  less  than  5  m's  model  (Noonkester,  1985).  Here,  the  size 
distribution  is  not  represented  by  a  series  of  lognormals  (as  in  Figure  3),  ‘here  is 
no  requirement  for  calculations  to  fit  NAM  predictions  at  10  meters,  and  the 
wavelengths  available  for  NOVAM  calculation  are  from  1  to  II  nm.  The  limita¬ 
tion  on  wavelengths  is  due  to  optical  calculations  using  a  simplification  of  .Mie 
scattering  parameters. 

3.  A  weak  convection  model  (Davidson  &  Fairall,  1986)  used  when  scattered  cumulus 
clouds  are  present,  cloud  tops  do  not  exceed  3  km,  and  a  well-mixed  layer  is  present 
below  the  cloud  base, 

4.  A  default  profile  (Gathman,  1989)  which  is  used  when  there  is  no  vertical  potential 
temperature  and  mixing  ratio  information. 

A  stratus  model  for  winds  greater  than  5  m;s  and  a  deep  convection  model  for  cloud  tops 
extending  above  3000  m  arc  not  yet  supported. 
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Figure  5.  Stylized  proFile  definitions  for  use  in  preamble:  (Gatiunan,  1989). 


NOVAM  sub-model  selection  is  based  on  the  input  parameters  describing  the  ver¬ 
tical  stratification  (thermal  stability,  presence  of  an  inversion,  and  the  inversion  height), 
cloud  cover,  cloud  type,  wind  speed,  and  the  requested  wavelength  for  extinction  calcu¬ 
lations. 

1.  The  default  profile  generator 

The  default  profile  generator  uses  the  surface  meteorological  measurements  to 
develop  profiles  of  temperature  and  relative  humidity,  generating  values  for  every  100 
meters  of  altitude.  The  default  profile  differs  from  a  radiosonde  derived  profile  by  having 
the  vertical  distribution  of  the  lognormal  components  described  by  an  exponential 
function  of  scale  height.  For  each  default  altitude  level,  the  optical  calculations  arc  based 
on  the  combined  size  distributions  of  the  aerosol  deduced  from  the  model  assumptions. 
The  index  of  refraction  and  size  arc  determined  by  the  response  (swelling)  of  hygroscopic 
aerosols  based  on  the  relative  humidity  calculated  from  the  relative  humidity  profile 
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Figure  6.  NOVAM  flow  diagram:  (Gathinan,  1989). 

generator  (Gathman,  78).  I'roni  these  parameters,  the  aerosol  extinction  at  each  level 
can  be  calculated  directly. 
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IV.  INFLUENCES  ON  AEROSOL  SIZE  DISTRIBUTION  PROFILES 

WITH  WEAK  CONVECTION 


A.  CONCENTRATION  GRADIENTS  IN  NOVAM 

The  basis  for  modeling  the  evolution  of  aerosol  concentration  is  the  continuity 

equation.  In  the  Boussinesq  approximation,  which  assumes  a  shallow  planetaiy 

dp 

boundarv'  layer,  an  incompressible  fluid  =  0).  and  negligible  horizontal  gradients  of 
density,  an  aerosol  component  concentration,  X,  can  be  written  in  its  simplest  one- 
dimensional  form  as 


+  IVjjX  +  IE  4^ 
Cl  vz 


cz  *  cz 


where  L'  is  the  mean  horizontal  wind. 

\V  is  the  mean  vertical  wind, 

is  the  particle  gravitational  settling  velocity, 

D  is  the  particle  molecular  transpo’".  coefTtcient, 

IVX  is  the  mean  vertical  flux  of  particles  due  to  turbulent  transport. 


(7) 


In  this  equation,  X  must  be  a  conservative  aerosol  variable,  e.g.  dr>'  size  concentration, 
that  is  unalTccted  by  variations  in  relative  humidity.  Terms  such  as  If'.V  are  unlikely  to 
be  known  from  first  principles  leading  to  applying  models  of  the  boundaiy  layer  to  sim¬ 
plify  the  process.  (Davidson  and  Fairall,  I9S6) 

The  .Albrecht  model  was  developed  specifically  for  a  trade-wind  region  (.Albrecht, 
1979)  and  has  been  incorporated  in  NOV.AM  to  describe  the  weak  convective  regime. 
The  three-layer  structure  of  this  regime,  described  by  .Albrecht,  is  shown  in  Figure  7. 
In  Figure  7,  C  represents  the  dry  aerosol  concentration  with  the  subscripts  r,  b,  c,  and 
i  representing  the  surface,  cloud  base,  the  in-cloud  layer,  and  the  cloud  top  respectively. 
The  layers  and  their  characteristics  arc  summarized  below,  Davidson  &  Fairall  (1986): 

1.  The  subcloud  layer  (0  <  z  <  r*)  where  is  the  cloud  base  (LCL).  This  layer  is 
characterized  by  linear  fiux  profiles  and  '  well-mixed"  mean  profiles.  The  depth  of 
this  layer  is  changed  by  subsidence  and  entrainment. 

2.  The  cumulus  cloud  layer  (Zj  <  z  <  zj.  characterized  by  parabolic  (lux  profiles  and 
linear  mean  profiles,  'fhe  slopes  of  the  mean  properties  are  determined  by  empirical 
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models  of  cloud  transport  and  entrainment  properties  as  well  as  the  boundar>' 
conditions  at  the  top  and  bottom  of  the  cloud  layer. 

3.  I'he  free  troposphere  (non-turbulent)  above  the  capping  "trade"  inversion.  The 
properties  of  this  layer  are  determined  by  advcction  and  are  usually  input  variables 
to  the  model. 


Note  that  the  mean  profiles  are  not  continuous  at  interfaces  2*  and  2,. 

Factors  that  influence  the  in-cloud  dry  aerosol  concentration  gradient,  y  (shown  in 
Figure  7),  are  given  by  the  following  balance  equation: 


_£l 

Di 


=  (V.F-  )y  -h 


(-.-2,)' 


(8) 


where  V  V  is  divergence, 

iv‘,  is  the  entrainment  rate  of  the  inversion, 
t  is  the  cloud  relaxation  time  in  days, 

E  is  the  cloud  wall  entrainment  variable  given  by: 


E  = 


(*/ 


2fthV-t-2A7;(l  -fF) 


eb 


(9) 


where  AT,  =  .5'K. 

r:^1.6, 

y,  =  the  entrainment  rate  of  0,  in  the  layer, 
is  the  jump  in  0,  at  the  cloud  base. 

Because  the  layer  below  the  cloud  base  is  considered  well-mixed,  the  dry  aerosol  con¬ 
centration  gradient  for  this  layer  is  normally  considered  to  be  0,  but  a  non-zero  gradient 
could  be  predicted  by  the  model  with  the  specified  parameters.  (Davidson  and  Fairall, 
I9S6) 

Some  assumptions  were  made  to  simplify  calculations  when  incorporating  the 
Albrecht  model  into  NOVAM.  The  Albrecht  model  is  dynamic,  predicting  the  temporal 
evolution  of  the  boundary  layer.  NOV.AM  seeks  to  characterize  the  average  or  typical 
vertical  structure,  therefore  all  time  derivatives  are  neglected  in  the  above  equations. 
.Additionally,  the  following  assumptions  are  used  to  simplify  calculations:  (Davidson  & 
Fairall.  1986) 

I.  In  the  mi.xcd  layer.  C  =  C 
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Figure  7.  Schematic  diagram  of  aerosol  concentration:  Dry  aerosol  concentration 
based  on  maritime  and  continental  sources  for  a  weak  cumulus  con¬ 
vection  boundary  layer  structure.  (Gathman,  I9S9) 

2.  In  the  cloud  layer,  C=  C;  +  y,(r  —  r*) 

3.  C,  is  given  by  a  standard  wind  speed  dependent  model 

4.  At  the  cloud  top,  assume: 

a.  C;  =  0  if  C  represents  a  locally  generated  sea-salt  mode 

b.  C*  is  given  by  climatology 

5.  Q  =  C;  =  C, 

These  assumptions  lead  to  the  following  equation  for  the  vertical  gradient  from 
cloud  base  to  cloud  top: 


(Ct-Q)n 
( I  +  /i)A* 


(10) 


where  n  is  a  dimensionless  entrainment  parameter  at  the  cloud  top  given  by: 


n  = 


4vv*^x  (t-V-n 
Az 


(11) 


In  this  form,  the  entrainment  rate  at  the  cloud  base  and  tne  cloud  wall  entrainment  are 
not  important.  These  simplifications  may  reduce  the  accuracy  in  the  vertical  represen¬ 
tation  of  dry  aerosol  concentration,  but  this  predicted  vertical  structure  is  a  large  first 
step  in  improving  the  exponential  model  currently  in  use.  With  mixing  and  entrainment 
controlled  by  the  potential  temperature  and  relative  humidity  profiles,  a  more  accurate 
representation  of  the  transport  properties  in  the  sub-cloud  and  cloud  layers  is  obtained. 

B,  RELATIVE  HUMIDITY  EFFECTS  IN  NOVAM 

Since  aerosol  particles  attenuate  energy  in  an  amount  proportional  to  their  cross 
sectional  area,  any  model  requiring  aerosol  extinction  calculations  must  take  into  ac¬ 
count  the  elfect  of  relative  humidity  on  aerosol  growth.  Vertical  changes  in  relative 
humidity  become  important  when  requiring  slant-path  aerosol  extinction  calculations  for 
the  weak  convective  region  where  relative  humidity  is  expected  to  increase  with  height 
up  to  the  cloud  layer,  then  decrease  gradually  to  top  of  the  cloud  layer.  Figure  8  shows 
the  effect  relative  humidity  has  on  particle  radius  (represented  by  the  ratio  of  the  aerosol 
particle  radius  at  the  indicated  relative  humidity  (r)  versus  the  radius  of  the  same  aerosol 
particle  at  0%  relative  humidity  (rj).  Fitzgerald  found  that  pure  NaCl  particles  experi¬ 
enced  a  sudden  increase  in  size  at  some  critical  value  of  relative  humidity  (between  70 
and  76%  relative  humidity  as  represented  by  the  hatched  area  in  Figure  8)  wiiereas  na¬ 
turally  occurring  maritime  aerosol  particles  show  continuous  smooth  growth  with 
changing  relative  humidity.  This  may  be  explained  by  the  gradual  deliquescence  of  sol¬ 
uble  material  which  is  present  as  a  mixture  of  different  salts. 

Based  on  the  information  shown  in  Figure  8,  Fitzgerald  developed  the  following 
relationship: 


riS) 


(12) 


where  y  represents  the  air-mass  characteristic  and  S  represents  the  saturation  ratio  (for 
the  purposes  of  this  discussion,  S  can  be  assumed  to  be  equal  to  RH  /  100). 

The  relative  humidity  dependence  in  both  NAM  and  NOVA.VI  appears  in  the  form 
of/  in  the  equation  for  size  distribution  (Equation  3).  As  can  be  seen  in  Figure  8,  when 
the  relative  humidity  increases,  deviation  of  the  measured  ratios  from  any  theoretical 
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Figure  8.  Growth  curves  for  marine  aerosols:  Theoretical  and  calculated  growth 
curves  for  NaCl  particles  and  natural  marine  aerosols  (Fitzgerald,  1978). 
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curve  which  is  based  on  dry  particle  size  becomes  greater.  Because  relative  humidity  is 
generally  above  50%  throughout  the  MABL,  a  reference  relative  humidity  of  80%  was 
selected  to  minimize  the  error  inherent  in  the  theoretical  aerosol  growth  calculation,  and 
this  reference  relative  humidity  was  used  as  a  basis  in  computing/.  An  equation  for/ 
is  generated  by  dividing  equation  12  by  the  ratio  of  r(.8)  to  r,  yielding  (Gerber,  1985): 


ft  = 


r{S) 

;(.8) 


a  I  -  s 

C8i(l  -  S) 


1 

3 


(13) 


where  C7,  C8  are  conditional  constants  based  on  the  aerosol  mode  component  (recall 
Figure  3  on  page  7)  as  displayed  in  Table  3.  To  incorporate  the  fourth  mode  (mode  0) 
into  Equation  3,  simply  change  the  summation  to  include  modes  0  through  3. 
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Table  3.  GROWTH  FACTOR  CONSTANTS  AND  VALIDITY  RANGE 


Aerosol  Model 
Composition 

Mode 

Component 

(i) 

C7, 

C8, 

Range  of 
Validity 

Sea  salt 

3 

1.97 

5.83 

RH  <  99.99% 

Sea  salt 

2 

1.83 

5.13 

RH  <  99.9% 

Water  Soluble 

1 

1.17 

1.87 

RH  <  99% 

Dust 

0 

/o^I.O 

all  values 

% 
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V.  THE  KEY90  EXPERIMENT 


A.  GOAL 

Most  observational  studies  within  the  development  of  NOVAM  have  been  accom¬ 
plished  off  the  California  coast.  Initial  experiments  involving  the  testing  of  NOVAM 
were  performed  in  the  same  environment  that  the  model  was  developed  in.  As  expected, 
the  model  verifies  reasonably  well  in  that  environment.  Since  NOVAM  is  intended  to 
be  a  global  maritime  model,  the  model  had  to  be  tested  in  other  geographical  scenarios. 

The  KEY90  experiment  was  the  first  attempt  to  test  N’OVAM's  performance  in  a 
tropic-like  environment  where  shallow  or  deep  convection  exists.  The  shaded  region  in 
Figure  1,  known  as  Marathon,  was  the  region  of  data  collection  during  KEY90. 
Insolation  heating  of  the  relatively  shallow  waters  between  Florida  and  Cuba  causes 
convection  to  be  the  major  meteorological  process  taking  place  near  the  air-sea  interface 
in  this  region.  Also,  at  .Marathon  data  could  be  obtained  by  both  boat  and  aircraft  while 
being  away  from  land  influences  and  major  continental  effects  on  the  data. 

B.  MEASUREMENTS 

Simultaneous  measurements  were  made  of  parameters  that  the  model  requires  as 
inputs  and  of  extinction  profiles  at  various  wavelengths.  The  extinction  profile  meas¬ 
urements  are  then  compared  to  the  results  obtained  by  running  N'OVA.M  at  wavelengths 
corresponding  to  those  measured.  The  extinction  profile  measurements  provided  the 
standard  to  judge  NOVAM  extinction  profile  predictions.  Although  this  procedure 
seems  fairly  straight  forward,  it  is  noted  that  it  is  difficult  to  provide  an  accurate  estimate 
of  either  tlie  measured  extinction  profile  or  of  the  quality  of  the  meteorological  data  used 
as  inputs  to  N'OVA.M  (e.g.  surface  IR  extinction  @  10.6  nm  and  a.m.p.). 

The  data  measurements  collected  during  KEY90  were  made  from  shore,  from  a  boat, 
and  from  two  aircraft.  Both  the  atmospheric  surface  layer  and  the  .MABL  structures 
were  probed  in  detail.  Figure  9  shows  the  instrumentation  involved  in  the  experiment. 
For  each  instrument,  an  asterisk  marks  the  days  of  operation.  Instruments  are  grouped 
according  to  the  information  provided  (i.e.  air  temperature,  sea  surface  temperature, 
relative  humidity,  wind  speed,  wind  direction,  atmospheric  pressure,  visibility,  vertical 
profiles,  Lidar  information,  and  other  data).  Table  4  gives  a  summary  of  the  time  and 
location  that  boat  and  aircraft  atmospheric  vertical  profiles  were  taken.  For  easy  refer- 
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Figure  9.  Suniiuary  of  instrumentation  used  during  KEY90;  (Gatliman,  1991) 


ence,  the  difTerence  of  time  (in  minutes)  and  in  location  (in  nautical  miles)  is  provided 
in  the  last  two  columns. 


22 


Table  4.  SOURCES  FOR  VERTICAL  PROFILES 


Aircrafl 

Radiosonde 

Day 

Time 

(UTC) 

Location 

Time 

(UTC) 

Location 

At 

(min) 

RQI 

09  Jul 

2120 

A 

2200 

24M3'iV80°33'lF 

40 

23 

10  Jul 

2100 

24°39'iV80“55'lF 

2200 

24°23LV  80°33TF 

60 

25 

12  Jul 

0030 

B 

2330' 

24’21LY80'33'IF 

60 

26 

12  Jul 

1930 

B 

1906 

24M5'A' 80*21' IF 

24 

39 

1945 

24°21'A' 80*33' IF 

15 

39 

13  Jul 

1330 

B 

1319 

24*21'A' 80*33' IF 

11 

26 

14  Jul 

1140 

B 

1105 

24*21'A' 80*33' IF 

35 

26 

16  Jul 

2100 

A 

2027 

24*13'A'80“31'1F 

33 

25 

2359 

B 

2324 

24*21'A' 80*32' IF 

35 

27 

17.1ul 

.  _  . 

1605 

B 

1617 

24*36'A' 80*53' IF 

12 

2 

18  Jul 

1240 

B 

1218 

24*30'A' 80*53' IF 

22 

5 

1545 

A 

1459 

24*30'A'8r00'lF 

46 

9 

19  Jul 

1020 

B 

1023 

24*30'A' 80*53' IF 

3 

5 

'  11  Jul 

A  =  24“22'/V80‘’55'lF 
B  =  24“35'iY80‘’55'lK 


It  should  be  noted  that  in  the  process  of  collecting  data,  areas  of  thunderstorm  ac¬ 
tivity  were  avoided.  This  may  have  caused  parameters  such  as  average  wind  to  be  lower 
than  what  may  have  actually  o  -curred. 

C.  ATMOSPHERIC  SYNOPTIC  SCALE 

The  major  synoptic  flow  patterns  w'hich  influenced  conditions  in  the  .Marathon  re¬ 
gion  from  9  to  19  July  are  shown  in  Figure  10.  For  the  period  of  9  -  10  July,  surface 
winds  at  Marathon  were  generally  from  the  east  to  south-east  with  speeds  of  approxi¬ 
mately  5  m/s.  This  was  a  result  of  the  influence  of  high  pressure  centered  off  the  east 
coast  of  the  U.S.  (Figure  10a).  The  high  pressure  system  continued  to  control  the 
Marathon  region  winds  on  the  11th  and  12th  of  July,  however,  short  wave  troughs  were 
analyzed  as  passing  through  Marathon  on  11  July. 

As  the  high  pressure  system  moved  to  the  north-east,  winds  at  Marathon  shifted 
from  east-southeast  to  a  more  southerly  direction  (Figure  10b).  Wind  speeds  were  be- 


F.^ure  10.  Synoptic  pattern  in  Marathon  area:  Surface  pressure  (mb)  and  700 
mb  heights  (m  *  10'‘  )  from  NMC  analysis  for  1200  UTC  is  provided 
for  the  following  days:  a)9  July  b)12  July. 
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Figure  1 1.  Synoptic  pattern  in  Marathon  area  (cont.):  Surface  pressure  (mb)  and 
700  mb  heights  (m  *  10-'  )  from  NMC  analysis  for  1200  UTC  is  pro¬ 
vided  for  the  following  days:  a)  1 4  July  b)  17  July. 
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tween  3  and  5  m's  during  this  period.  By  13  July,  a  low  pressure  system  which  was  ap¬ 
proaching  the  eastern  U.S.  had  an  influence  on  winds  at  Marathon.  Wind  flow  remained 
southerly,  but  wind  speed  decreased  slightly  due  to  weaker  pressure  gradients  over  the 
area.  On  14  July,  the  low  pressure  center  had  deepened  and  moved  to  the  east 
(Figure  11a).  Surface  wind  flow  over  Marathon  remained  southerly  with  w'ind  speeds 
increasing  to  6  -  8  m/s.  The  cold  front  associated  with  the  low  pressure  system  was 
moving  very  slowly  to  the  east.  By  late  UTC  on  this  day,  it  was  evident  that  the  frontal 
system  was  losing  its  upper  level  support,  and  the  high  pressure  evident  at  700  mb  near 
the  Florida  Keys  in  the  figure  continued  to  build. 

From  15  to  19  July,  high  pressure  once  again  dominated  the  flow  pattern  in  the 
Marathon  region  (Figure  1  lb).  As  the  high  pressure  center  moved  to  the  north-east,  the 
pressure  gradient  at  all  levels  weakened.  Surface  winds  were  from  the  east  to  south-east 
throughout  this  period  with  wind  speed  decreasing  to  approximately  5  m/s  by  the  19th. 
A  series  of  short-wave  troughs  were  analyzed  as  passing  th:  uugh  the  Marathon  region 
on  17  July. 

Although  all  observations  during  KEY90  fell  within  tiie  realm  of  weak  convection 
defined  by  NOVA.M,  there  were  examples  of  deep  convection  in  the  form  of 
thunderstorms.  Thunderstorm  activity  was  observed  in  the  Marathon  region  on  14,  16, 
17,  and  19  July. 


% 
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VI.  NOVAM  USAGE,  PROCEDURE  AND  RESULTS 


A.  PROCEDURE 

1.  Surface  File  Generation 

Each  radiosonde  and  flight  spiral  profile  requires  corresponding  surface  file. 
Figure  9  shows  that  there  were  many  choices  for  most  surface  file  parameters.  Because 
of  consistent  performance  during  the  entire  experiment,  WeatherPak  measurements  were 
used  for  surface  layer  air  temperatures,  wind  speeds,  and  relative  humidities.  When  this 
information  was  not  available  on  10  and  1 1  July,  sling  psychrometer  measurements  were 
used  for  air  temperature  and  relative  humidity  while  a  hand  held  device  was  used  for 
wind  speed.  Bucket  thermometer  measurements  of  sea  surface  temperature  were  used. 
NRL's  radon  counter  was  used  to  establish  the  air  mass  parameter.  The  NOSC  aircraft 
aerosol  spectra  was  used  to  obtain  the  10.6  fxm  aerosol  extinction  at  the  surface.  The 
altitude  of  this  reading  was  usually  30  meters  above  the  surface.  The  optical  visibility 
value  was  calculated  by  the  model,  based  on  the  above  parameters. 

2.  Radiosonde  Preamble  Generation 

The  general  shapes  of  the  potential  temperature  and  mixing  ratio  (for  flight)  or 
specific  humidity  (for  radiosonde)  profiles  were  examined  together  when  generating  the 
radiosonde  preamble.  The  first  step  was  to  establish  heights  of  the  first  and  second  inv¬ 
ersions  (2*  and  2,  in  Figure  5).  In  general,  a  well-mixed  layer  was  easy  to  establish  below 
the  first  inversion,  where  potential  temperature  and  ntixing  ratio  remained  constant.  The 
more  difficult  task  was  determining  a  height  for  the  second  inversion.  In  most  cases, 
there  was  no  clear  boundary  (jump)  to  define  the  top  of  the  cloud  layer,  especially  in  the 
potential  temperature  profiles.  As  a  result,  relative  humidity  profiles  and  observations 
of  slope  changes  in  the  potential  temperature  and  mixing  ratio  profiles  were  considered 
subjectively  in  making  the  determination. 

3.  NOVAM  Runs 

NOVAM  runs  for  both  the  aircraft  and  radiosonde  vertical  profiles  used  the 
same  surface  file,  so  only  one  surface  file  was  created  for  each  row  presented  in 
Table  4.  Information  in  the  surface  file  is  displayed  above  the  radiosonde  vertical  pro¬ 
file  panels  in  Figure  14  through  Figure  26. 

The  jumps  occurring  between  straight  line  segments  of  the  radiosonde  and  air¬ 
craft  profiles  are  important  run  parameters  The  predicted  amount  of  entrainment  oc- 
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curring  at  the  cloud  top  (second  inversion)  is  dependent  on  the  difference  (jumps)  in  the 
potential  temperatures  and  mixing  ratios  at  the  level.  The  smaller  the  difference,  the 
more  entrainment  is  occurring  between  the  cloud  layer  and  the  layer  above  the  second 
inversion.  The  wavelength  of  interest  for  each  model  run  was  an  input  parameter  since 
aerosol  extinction  varies  with  wavelength.  All  the  results  presented  in  this  thesis  are  for 
a  wavelength  of  3.5  urn,  even  though  aerosol  extinction  profiles  w'ere  obtained  by  NOSC 
aircraft  at  .53  urn,  1.06  urn,  3.5  fini,  and  10.6  nm. 

A  NOVAM  output  file  was  generated  for  each  vertical  profile  presented  in  Table 
4.  A  separate  run  was  also  made  to  generate  a  limiting  profile.  For  each  surface  file,  the 
XOVA.M  default  vertical  profile  was  used  with  the  relative  humidity  at  each  level  set  at 
98%.  When  aerosol  extinction  measurements  in  the  low'er  portion  of  the  MABL  exceed 
this  value,  precipitation  is  believed  to  be  a  possible  cause,  and  the  assumption  of 
NOVAM  that  no  precipitation  exists  is  violated.  The  assumption  of  the  default  profile 
generator  that  each  mode  of  the  size  distribution  decreases  exponentially  with  height 
will  cause  a  difference  in  NOVAM  versus  limiting  profile  size  distributions,  increasing 
as  the  height  of  the  well-mi,xed  layer  increases. 

The  importance  of  the  relative  humidity  to  aerosol  e.xtinction  can  be  shown  by 
comparing  NOVAM  profiles  to  limiting  profiles.  At  a  wavelength  of  3.5  fjim,  an  increase 
in  relative  humidity  from  the  80%  typically  observed  at  the  surface  to  a  98%  relative 
humidity  (as  used  in  the  limiting  profiles)  results  in  an  approximate  order  of  magnitude 
increase  in  aerosol  extinction  prediction  (comparison  of  aerosol  extinction  at  surface 
between  NOVA.M  predictions  based  on  vertical  profiles  and  NOVAM  limiting  profiles 
in  Figure  14  through  Figure  26).  At  the  times  of  the  vertical  profiles,  relative  humidities 
above  90%  were  rarely  seen,  therefore  a  comparison  could  not  be  made  of  how'  well 
NOVA.M  actually  performed  at  the  higher  relative  humidities. 

B.  RESULTS 

Simple  comparisons  between  the  flight  and  radiosonde  potential  temperature,  mix* 
ing  ratio,  and  aerosol  extinction  profiles  reveal  extensive  variability  in  space  and  time. 
Results  from  19  July  (Figure  26)  illustrate  this  point  the  best.  From  Table  4,  the  ap¬ 
proximate  horizontal  distance  between  vertical  profiles  was  five  nautical  miles,  while  the 
time  difference  was  negligible.  Analyses  to  determine  the  cloud-layer  boundaries  (defined 
by  the  first  and  second  inversion)  yielded  a  150  meter  difference  in  height  of  the  cloud 
base  and  a  700  meter  difference  in  the  height  that  the  second  inversion  occurred.  The 
fact  that  there  was  a  substantial  spatial  variability  is  further  corroborated  by  the  Lidar 
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extinction  profile  for  14  July  shown  in  Figure  27,  which  clearly  shows  the  spatial  vari¬ 
ation  in  atmospheric  structure.  As  discussed  above,  slant-path  calculations  of 
transmittance  currently  assume  aerosol  extinction  to  be  horizontally  uniform.  This 
condition  does  not  appear  to  exist  in  the  KEY90  region.  .A  predicted  aerosol  extinction 
profile  generated  by  NOVAM  based  on  one  vertical  profile  will  most  likely  not  be  rep- 
*■  resentative  of  the  average  vertical  structure  over  the  path  of  interest.  In  order  for 

NOVAM  to  provide  the  required  predictions,  either  the  mean  cloud  base  and  cloud  top 
*  levels  over  the  path  need  to  be  established,  or  a  technique  needs  to  be  developed  to 

minimize  the  importance  of  the  analyzed  inversion  levels. 

All  NOSC  aircraft  aerosol  extinction  observations  between  9  -  19  July  were  plotted 
simultaneously  (Figure  12)  to  e.xamine  how  well  NOVA.M  predictions  characterized 
aerosol  extinction.  A  subjectively  determined  mean  profile  was  generated  from  the  plots 
in  Figure  12.  and  will  be  referred  to  as  the  flight  composite.  To  be  considered  successful, 
separate  NOVAM  predictions  should  compare  favorably  to  the  flight  composite  under 
most  circumstances.  Due  to  space  and  time  differences  between  the  flight  profiles  and 
the  radiosonde  profiles,  only  NOVAM  results  from  flight  profiles  were  used  n  t/ie  com¬ 
parisons  with  the  observed  aerosol  extinction  profiles.  When  considering  only  the  top 
row  of  Figure  14  through  Figure  26,  the  following  features  can  be  seen: 

1.  NOVA.M  aerosol  extinction  predictions  above  the  second  inversion  were  extremely 
poor. 

2.  Below  the  second  inversion,  N’OV^AM  predictions  compared  favorably  with  the 
flight  composite  except  in  the  cases  where  aerosol  extinction  dropped  off  noticeably 
(an  order  of  magnitude  difference  between  the  surface  and  cloud  base  observations) 
below  the  cloud  base  (see  top  panel  Figure  17c.  Figure  20c,  Figure  24c,  and  F'g- 
ure  25c). 

The  assumption  that  aerosol  extinction  above  the  second  inversion  is  essentia'  y 
zero  when  the  air  mass  parameter  is  calculated  to  be  approximately  equal  to  one  d»  s  s 
not  appear  to  be  a  good  one  for  this  region.  Aerosol  extinction  was  determined  to  (ve 
significant  if  it  is  within  one  order  of  magnitude  of  the  aerosol  extinction  that  NOVM  >  ' 
'  calculates  at  the  surface  for  the  wavelength  of  interest.  Using  this  as  a  criteria,  the  'i. 

13,  14,  17,  and  18  July  NOVA.M  flight  profiles  all  showed  significant  aerosol  cxtinct;cn 
»  contributions  above  the  second  inversion. 

One  possible  explanation  for  the  poor  performance  of  NOV.A.M  above  the  second 
inversion  is  that  the  weak  or  shallow  convection  assumptions  arc  not  cflcctivc  here  .i.e. 
this  may  be  a  case  of  deep  convection.  Thunderstorm  activity  and  rain  were  frequently 
observed  in  .Marathon,  and  multiple  cloud  layers  were  not  uncommon.  The  hypothesis 
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Figure  12.  Flight  extinction  profiles  nitli  flight  composite: 


is  that  relatively  small  regions  of  deep  convective  activity  result  in  an  increased  pro¬ 
duction  of  aerosol  at  the  surface  and  transports  aerosol  to  higher  levels.  As  these  deep 
convective  cells  dissipate,  large  amounts  of  aerosol  are  still  being  transported  horizon¬ 
tally  throughout  the  region,  which  would  account  for  the  aerosol  present  above  the  • 

second  inversion  in  many  cases. 

Another  possible  c.xplanation  may  be  that  the  atmosphere  above  the  second  inver¬ 
sion  conunonly  contains  a  significant  number  of  small  radii  aerosol  particles.  These 
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Figure  13.  NOVAM  prediction  profiles  with  NOV  AM  composite: 


particles  may  originate  either  from  continental  sources  or  from  the  ocean  surface  in  the 
form  of  DMS.  These  D.MS  particles  would  be  '•arried  to  the  upper  atmosphere  by 
convective  cells  and  transported  by  large  scale  circulation.  Particle  processes  associated 
with  DMS  origins  arc  not  well  understood.  {Durkee  ct  al.,  1990) 

If  the  predicted  aerosol  e.xtinction$  above  the  second  inversion  arc  ignored,  and  a 
simultaneous  plot  of  NOVA.M  predictions  based  on  flight  vertical  profiles  is  made,  a 
NOV.A.M  composite  can  be  obtained  (Figure  13)  in  the  same  manner  as  the  flight 
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composite  described  above.  The  flight  and  NOVAM  composites  turned  out  to  be  nearly 
identical  (Figure  14c).  There  was  not  much  variation  in  wind  speed  at  Marathon  and 
no  frontal  systems  passed  through  during  the  KEY90  experiment,  therefore  it  is  likely 
that  aerosol  generation  and  transport  mechanisms  remaiined  fairly  constant  during 
KHY90.  This  would  explain  why  composite  mean  profiles  performed  as  well  as  NOVAM 
predictions  below  the  second  inversion.  It  is  possible  that  mean  profiles,  such  as  the 
composites  generated  above,  can  characterize  the  typical  aerosol  extinction  (at  a  3.5 
urn  wavelength)  for  the  Marathon  region  during  July.  This  approach  leads  to  the  region 
being  described  empirically,  therefore  the  predictions  are  specific  to  one  region  during 
one  season.  The  required  data  base  for  describing  all  ocean  regions  for  each  season 
would  be  enormous  and  the  time  and  expense  involved  would  be  substantial. 

The  fact  that  a  simple  exponential  curve  developed  from  composite  mean  profiles 
performed  much  better  than  NOVAM  above  the  evaluated  second  inversion  indicates  a 
need  for  farther  research  on  aerosol  size  distributions  in  the  upper  atmosphere.  If  the 
structure  observed  in  KEY90  is  in  fact  due  to  deep  convective  processes,  a  revision 
would  have  to  be  made  to  the  sub-model  selection  process  in  NOVAM.  Also,  instead 
of  the  weak-convective  assumption  that  the  mode  2  and  3  components  drop  off  to  0 
above  the  second  inversion,  an  exponential  function  could  be  used  to  describe  the  aer¬ 
osol  size  distributions  for  these  modes. 

1 

It  is  believed  that  the  observations  where  aerosol  extinction  decreased  'noticeably 
below  the  first  inversion  are  attributed  to  the  combination  of  increased  aerosol  pro¬ 
duction  at  a  time  prior  to  the  observation,  possibly  caused  by  thunderstorm  activity,  and 
scavenging  in  the  lower  layer  due  to  rain,  creating  a  deficit  of  aerosol  below  the  cloud 
layer.  If  future  experiments  reveal  that  rain  scavenging  accounts  for  reduced  aerosol  ex¬ 
tinction  below  the  cloud  layer,  the  present  weather  surface  input  file  parameter  is  cur¬ 
rently  the  only  means  that  this  information  can  be  relayed  to  NOVAM.  Because  the 
aerosol  extinction  predicted  at  the  surface  is  tied  to  surface  observations  of  aerosol  ex¬ 
tinction,  the  task  of  modeling  this  phenomenon  would  still  remain. 
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Figure  14.  9  .luly  results;  The  top  panels  are  based  on  flight  vertical  profiles 

observed  09  July,  2120  UTC.  Bottom  panels  are  based  on  radiosonde 
profiles  for  09  July,  2200  UTC.  Panels  labeled  a)  represent  temperature 
and  mixing  ratio  (flight)  or  specific  humidity  (sonde)  versus  altitude, 
b)  plots  relative  humidity  versus  altitude,  and  c)  shows  aerosol  ex¬ 
tinction  at  3.5  urn  versus  altitude.  Refer  to  text  for  a  description  of  the 
profiles  appearing  in  panel  c). 
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MIXING  RATO  (9/kq) 


Figure  15.  10  July  results;  The  top  panels  are  based  on  flight  vertical  profiles 

observed  10  July,  2100  UTC.  Bottom  panels  are  based  on  radiosonde  t 

profiles  for  10  July,  2200  UTC.  Panel  description  is  same  as  in 
Figure  14. 
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Figure  16.  12  July  results  for  0030  UTC:  The  top  panels  are  based  on  flight 

vertical  profiles  observed  12  July,  0030  UTC.  Bottom  panels  are  based 
on  radiosonde  profiles  for  1 1  July,  2330  UTC.  Panel  description  is 
same  as  in  Figure  14. 
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Figure  17.  12  July  results  for  1930  UTC  (a):  Tlie  top  panels  are  based  on  flight 


vertical  profiles  observed  12  July,  1930  UTC.  Bottom  panels  are  based 


on  radiosonde  profiles  for  12  Jul,  1906  UTC.  Panel  description  is  same 


as  in  Figure  14. 
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Figure  18.  12  July  results  for  1930  UTC  (b):  The  top  panels  are  based  on  flight 

vertical  profiles  observed  12  July,  1930  UTC.  Bottom  panels  are  based 
on  radiosonde  profiles  for  12  Jul,  1945  UTC.  Panel  description  is  same 
as  in  Figure  14. 
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Figure  19.  13  July  results:  The  top  panels  are  based  on  flight  vertical  profiles 

observed  13  July,  1330  UTC.  Bottom  panels  are  based  on  radiosonde  « 

profiles  for  13  July,  1319  UTC.  Panel  description  is  same  as  in 
Figure  14. 
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Figure  20.  14  July  results;  The  top  panels  are  based  on  flight  vertical  profiles 


observed  14  July,  1140  UTC.  Bottom  panels  are  based  on  radiosonde 


profiles  for  13  July,  1105  UTC.  Panel  description  is  same  as  in 


Figure  14. 


39 


Figure  21.  16  July  results  for  2100  UTC:  The  top  panels  arc  based  on  night 

vertical  profiles  observed  16  July,  2100  UTC.  Bottom  panels  are  based 
on  radiosonde  profiles  for  16  July,  2027  UTC.  Panel  description  is 
same  as  in  Figure  14. 
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Figure  22.  16  July  results  for  2359  UTC:  The  top  panels  arc  based  on  flight 

vertical  profiles  observed  16  July,  2359  UTC.  Bottom  panels  arc  based 
on  radiosonde  profiles  for  16  July,  2324  UTC.  Panel  description  is 
same  as  in  Figure  14. 
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Figure  23.  17  July  results  for  1605  UTC:  The  top  panels  arc  based  on  flight 

vertical  profiles  observed  17  July,  1605  UTC.  Bottom  panels  •_  ’.'a.sed 
on  radiosonde  profiles  for  17  July,  1617  UTC.  Panel  de.scnption  is 
same  as  in  Figure  14. 
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Figure  24.  18  July  results  for  1240  UTC:  The  top  panels  arc  based  on  flight 

vertical  profiles  observed  18  July,  1240  UTC.  bottom  panels  are  based 
on  radiosonde  profiles  for  18  July,  1218  UTC.  Panel  description  is 
same  as  in  Figure  14. 
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Figure  25.  18  July  results  for  1545  UTC;  The  top  panels  arc  based  on  flight 

vertical  profiles  observed  18  July,  1545  UTC.  Bottom  panels  arc  based  4 

on  radiosonde  profiles  for  18  July,  1459  UTC.  Panel  description  is 
same  as  in  Figure  14. 
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Figure  26.  19  July  results:  The  lop  panels  are  based  on  flight  vertical  profiles 

observed  19  July,  1020  U  TC.  Bottom  panels  are  based  on  radiosonde 
profiles  for  19  July,  1023  UTC  Panel  description  is  same  as  in 
Figure  14. 
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Figure  27.  Logaritlunic.il',y  ajnplifief.  NRL  Lidar  extinction  profiles:  Coded  in 
false  color,  horizontal  axis  is  distance  coverec  jy  NRL  aircraft  from 
1107  -  1126  LMC  recorded  on  14  July  1990.  Vertical  axis  is  altitude 
above  sea  level.  Colors  represent  aerosol  extinction  in  knr'  as  given  by 
tlie  bar  code  on  the  right  hand  side.  (Hooper,  1991) 


46 


VII.  SUMMARY,  CONCLUSION  AND  RECOMMENDATIONS 

An  experiment  was  conducted  in  the  Florida  Keys  from  2-19  July  1990  to  test  the 

I 

^  performance  of  the  Naval  Oceanic  Vertical  Aerosol  Model  (NOVAM)  in  a  tropical, 

weak-convective  regime.  Meteorological  data  collected  by  aircraft  and  boat  was  used  to 
V  generate  the  surface  and  vertical  profile  information  files  required  by  NOVAM.  Using 

this  information,  NOVAM  predicts  the  aerosol  extinction  {krrr^)  for  a  vertical  cross- 
section  of  the  atmosphere.  When  NOVAM  is  incorporated  into  the  LOWTRAN  model, 
which  also  accounts  for  molecular  extinction,  slant-path  calculations  for  transmittance 
may  be  performed. 

Aircraft-observed  aerosol  extinction  profiles  were  obtained  in  order  to  verify 
NOVAM  predictions  of  aerosol  extinction.  Comparisons  between  observed  and 
NOVAM  aerosol  extinction  profiles  revealed  major  deviations  above  the  cloud  top. 
From  the  surface  to  the  top  of  the  cloud  layer,  NOVAM  generally  did  an  excellent  job 
in  predicting  profile  shape,  with  the  magnitude  of  aerosol  extinction  tied  to  the  ex¬ 
tinction  matched  at  the  surface.  In  a  few  cases,  observed  extinction  increased  more  rap¬ 
idly  than  NOVAM  predicted  extinction  from  the  surface  to  the  base  of  the  cloud  layer. 
This  is  attributed  to  rain  scavenging  associated  with  thunderstorm  activity  in  the  area. 
Comparison  between  difierent  aerosol  extinction  profiles  over  time  and  space  revealed 
much  variation,  which  Lidar  profiles  of  atmospheric  structure  verified.  Due  to  these 
variations,  many  vertical  soundings  of  the  atmosphere  would  have  to  be  obtained  in  or¬ 
der  to  characterize  a  mean  atmosphere  for  the  region  of  interest.  Thunderstorm  activity, 
multiple  cloud-layers,  and  the  spatial  variation  in  the  atmospheric  structure  have  led  to 
a  hypothesis  that  deep-convection  was  responsible  for  the  major  differences  between 
observed  and  predicted  aerosol  extinction  profiles.  If  this  is  the  case,  a  simple  modifica- 
tion  to  the  w'eak-convective  model  may  be  made  to  apply  to  a  deep-convection  model. 

The  results  described  above,  clearly  show  that  NOVAM  is  not  yet  ready  to  be  used 
operationally  in  regions  where  convection  is  the  dominant  process.  The  following  are 
recommendations  for  consideration  in  the  improvement  and  development  of  NOVAM: 

1.  In  order  to  eliminate  the  bias  in  wind  speed  measurements  caused  by  avoiding  areas 
of  storm  activity,  wind  measurements  from  buoys  placed  strategically  in  thef-  area 
of  operations  would  be  invaluable  in  generating  the  24  hour  average  winds.  '' 
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2.  Conduct  experiments  in  the  MARATHON  area  and  other  areas  where  weak 
convective  conditions  are  expected.  The  data  and  results  from  these  experiments 
could  then  be  compared  with  the  data  and  results  from  the  KEY90  experiment. 

3.  Continue  to  explore  the  possibility  that  some  cases  that  are  currently  classified  as 
weak  convection  may  in  fact  be  examples  of  deep  convection.  If  this  hypothesis  is 
correct,  :he  parameters  used  to  differentiate  between  shallow  and  deep  convection 
need  to  be  modified. 

4.  Explore  the  possibility  of  using  satellite  information  to  improve  model  predictions. 
One  example  would  be  to  use  the  techniques  established  by  Durkee,  et  al.  (1986; 
1990)  to  determine  variations  in  optical  depth  (total  extinction)  over  a  horizontal 
path  which  could  then  be  incorporated  into  NOV^A.M. 

5.  A  comparison  needs  to  be  performed  between  what  the  current  model  is  producing 
(LOWTRAN.NAM)  against  what  NOV.AM,  incorporated  in  LOWTRAN,  predicts 
(transmittance  predictions  versus  measured  would  be  extremely  useful). 
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